Absence of Phase Stiffness in the Quantum Rotor Phase Glass 
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We analyze here the consequence of local rotational-symmetry breaking in the quantum spin (or 
phase) glass state of the quantum random rotor model. By coupling the spin glass order parame- 
ter directly to a vector potential, we are able to compute whether the system is resilient (that is, 
possesses a phase stiffness) to a uniform rotation in the presence of random anisotropy. We show 
explicitly that the 0(2) vector spin glass has no electromagnetic response indicative of a supercon- 
ductor at mean-field and beyond, suggesting the absence of phase stiffness. This result confirms our 
earlier finding (PRL, 89, 27001 (2002)) that the phase glass is metallic, due to the main contribution 
to the conductivity arising from fluctuations of the superconducting order parameter. In addition, 
our finding that the spin stiffness vanishes in the quantum rotor glass is consistent with the absence 
of a transverse stiffness in the Heisenberg spin glass found by Feigelman and Tsvelik (Sov. Phys. 
g ; JETP, 50, 1222 (1979). 

o 

1 I. INTRODUCTION 

t/3 \ Spin glasses are characterized by the freezing of local spins along random non-collinear directions. Because each 
-4-J , spin points in a preferred direction, locally spin rotational symmetry is broken. Nonetheless, globally rotational 
symmetry is preserved because spin glasses have no net magnetization. We consider here the 0(2) quantum rotor 
model where the exchange interactions are random. As this model is isotropic in rotor space, a global rotation of all of 
the rotors is an exact symmetry, even in the glass phase. Nonetheless, in the glass state, a global rotation of all of the 
spins around any axis generates a new state which is distinguishable from the original unrotated state. Because such 
uniform rotations are generated by the group SO{2) , the spin glass state breaks 5*0(2) symmetry. All such states are 
energetically degenerate as a result of the inherent isotropy in rotor space. As a result of the broken SO (2) symmetry, 
, O, . it is reasonable to expect that a massless bosonic mode should exist. 

In the strict sense, a physical system possesses a non-zero phase rigidity if upon a uniform rotation of the phase, 
the free energy increase is of the form, 
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I/"") , where p s is the spin or superfluid stiffness and 8 is the collective phase variable. Consequently a spin-wave mode with 
■ a dispersion lu = ±ck would be an experimental signature of a spin stiffness consistent with Eq. Experimentally, 

C*~) I however, no such mode has ever been found in either neutron scattering or thermal measurementsi*Siii on spin glasses. 
This failure might be attributed to that fact that over-damped modes and/or low energy excitations conspire to make 
p s undetectable. Theoretically, in the phenomenological hydrodynamic account, Halperin and Saslow^ assumed that 

Sp s 7^ 0. They did caution the reader that the existence of a stiffness in a spin glass is subtle and, in all likelihood, 
doubtful as a result of the preponderance of experimental evidence for a large density of low-energy excitations that 
could over-damp the spin-wave mode. This conclusion is supported by extensive numerical simulations by Walker 
and Walstedt^ who found no evidence for the characteristic ui 2 vanishing of the low-energy modes. Two microscopic 
calculations of the spin stiffness exist. Feigel'man and Tsvelik^ developed a real-time diagrammatic technique for the 
Heisenberg spin glass and showed explicitly that the spin stiffness vanishes. This result is particularly robust because 
^ ' it follows from a simple permutation symmetry of the spin correlators^. Within the replica formalism of a Heisenberg 
spin glass Kotliar, Sompolinsky, and Zippelius 8,9 formulated a mean-field description of the single-valley stiffness 
• , constant. This limit is relevant at sufficiently short times that the spin glass remains trapped in a single configuration. 

In this limit, the stiffness constant is non-zero^. However, in the full statistical mechanical treatment of the problem 
in which hopping among the myriads of valleys in the energy landscape of a spin glass are allowed, the stiffness 
vanishes&Siifl. This result implies that the spin stiffness is a transient effect approaching zero in the equilibrium or 
long-time limit. In this limit, a new massless mode dispersing as k A emerges which leads to the vanishing of the spin 
stiffness, as in the real-time formalism^. Hence, there is a consilience between the replica and real-time formalisms 
that the stiffness constant vanishes in the Heisenberg spin glass. 

For quantum spin glasses, no calculation of the stiffness exists. Nontheless, we expect the same physics to be 
valid. Namely, as long as the system can relax and hop among all of the configurations of the spin glass, the stiffness 
should vanish. For example, in quantum spin glasses, quantum tunneling among the various local minima in the 
spin glass landscape is permitted, thereby leading to a vanishing of the stiffness. This problem is particularly current 
because we have recently proposed-* 1 that the bosonic excitations arising from fluctuations of the superconducting order 
parameter in the glassy phase, lead to a metallic conductivity at zero temperature. In the Gaussian approximation, this 
conductivity o^os diverges as 1/ m upon approaching the superconducting phase (m is the inverse correlation length of 
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the superconducting fluctuations). A free energy density of the form of Eq. Q, however, leads to a superconducting 
response. Hence, should the phase glass itself have a well-defined stiffness, then the bosonic conductivity, though 
intriguing, would be irrelevant as it would be dwarfed by the infinite conductivity arising from the excitations related 
to the glassy order parameter. We show here explicitly that this is not the case, at least at the mean-field level. 
Rather than attempting to calculate the phase stiffness from the free energy, we consider the linear response regime 
and couple the spin glass order parameter to the appropriate vector potential. Second, we compute the role of replica 
symmetry breaking (RSB) on the bosonic contribution to the conductivity. We show that weak RSB does not affect 
the metallic character of the conductivity as T — > 0. Consequently, the Bose metallic phase found earlier 11 is robust 
and constitutes the only known example of a metallic phase in 2D in the presence of disorder. 

II. PHASE STIFFNESS 

The starting point for our analysis is the 0(2) quantum rotor model, 

H = - E cJ2(w) -^JijcoBiOi-Oj-Aii), (2) 

where A y - = (e*/%) f. A • dl (e* = 2e). The Josephson couplings are assumed to be random and governed by a 
distribution 



p(Jy) = vib exp 



(Jjj ~ Jo) 

2J 2 



(3) 



with non-zero mean, Jo and J the variance. When the distribution has a non-zero mean, three phases are possible: 1) 
disordered paramagnet, 2) quantum phase glass, and 3) superconductor. Because the existence of the spin stiffness in 
the spin glass can be answered with the simpler model with zero mean (Jo = 0), we utilize this model at the outset. 
For a random system, the technique for treating disorder is now standard: 1) replicate the partition function, 2) 
perform the average over disorder and 3) introduce the appropriate fields to decouple the interacting terms that arise. 
As the corresponding action has been detailed previousljiiS^, we will provide additional steps that are necessary to 
determine how the electromagnetic gauge couples to the spin glass order parameter. We write the replicated partition 
function as 

= JVOiDJae- 8 (4) 

where the Euclidean action is given by 



I i,a 



where the superscript a represents the replica index. For J = 0, the integration over Jy in Eq. (0J results in the 
effective action, 
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+ yEE [ [ dTdT '\ E ^V^[et{r)~ae\{T')-{e a J {T)-ae){r'))-{A l3 {T)-aA l3 {T%} 

a,b (ij) •'° ^° a=+l,-l 
+C.C. (6) 

with a = +1,-1. As a result of the sum over a, we see that the vector potential enters both symmetrically and 
anti-symmetrically. To simplify the notation, we introduce the two-component vector 

S a (r) = (cos 6» a (r), sin <9 a (r)) (7) 

and the corresponding auxiliary field, 



Q;l(r,r') = (S;(r)Sl(r')) 



(8) 
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which will be used in decoupling the action and ultimately determines the Edwards-Anderson order parameter for 
the quantum spin glass transition. The remaining steps involve performing the cumulant expansion and taking the 
continuum limit. The final action can be separated into the local and gradient parts: 



iS'eff — Sloe + Sjr 



where the local part 



. Tl) T2)\ T1 = T2 = T 



-j J <iTl<iT 2 <iT a ^g; , 1 ,(x,T I ,T 2 )Qt=( X ,T 2 ,T a )g~( X ,T a ,T 1 ) 
a.b.c 

+ \ J ^^[ U Q-(x,r,r)g-(x,T ) T)+ U Q-(x ) r,T)Q-(x,r,r)] 
-| J d d x J dndrzYsiQ^n^)] 4 . 



(9) 
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is identical to that derived previously by Read, Sachdev and Ye** and the gradient part 



si- 
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V - ^-A(x, n ) + ^-A(x, r 2 ) ) Q^ 6 (x. r, . t-o ) 



V - ^-A(x, ri) - ^-A(x, r 2 ) j 0:; f 'i r, . r , i 
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in which the vector potential couples both symmetrically and asymmetrically to combinations of the Q— matrices of 
the same parity. Using the fact that Q± (t\, t 2 ) ~ (exp [i(9f (t) ± 6* b (t'))] ), the parity combinations of the Q— matrices 
are defined as follows: 



Q5?(x,t 1; t 2 ) = i [Q? 1 6 (x,t 1 ,t 2 )tQ^(x,ti,t 2 )] +^ [Q? 2 b (x,r 1 ,T 2 )±Q^(x,T 1 ,r 2 )] 



(12) 



It is evident that the vector potential enters in a non-time translationally invariant manner. This is a direct conse- 
quence of the fact that the Q— matrices themselves are a function of two independent times, not simply the difference 

Of Tl ~ T 2 . 

To calculate the conductivity, we need to focus entirely on the gradient part of the action as this is the only part 
that couples to the vector potential. The standard Kubo formula for the spin-glass contribution to the longitudinal 
conductivity takes the form, 



a(icj n ) = -— lim - [ d d (x - x') f d(r - r')—— 



S 2 Z n 



(x, t)SA x ( X ',t') 



iWn(T—T') 



(13) 



where we have chosen to orient the vector potential along the x-axis. A bit lengthy variational procedure leads to the 
following result: 



cr(iu) r , 



(e*) 2 1 



-£ f d(r-r')e^-^U f dr 2 ((|Q° b (x, r, r 2 )| 2 ) + (\Qf (x, r, r 2 )| 2 » S(r - r') 

U a.b J ° { h 



+4 «|Q?(x, r, r')| 2 ) - (|Q^ fc (x, r, r')| 2 » - / d d ( X - x')( J,(x, r) J x (x', /)) 



where the current J x (x, r) is defined as 



J(x,r) = l ^-Y, E / dr i Q^(x,r,r 1 )V(Q^(x,r,r 1 ))*-c.c 



ah a=+,— 



(14) 



(15) 



In deriving this expression for the current, we considered the relations <5" h (x, r 2 , ti) = Q" h (x, ti, t 2 ) and 
Q ab (x, t 2 , ri) = (Q_ (x, Ti,T2))*, that follow from the definition, Eq. I|12[l. To evaluate the correlation functions 
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in Eq. I|14|) . we need to use the Fourier components of the Q- fields: 

Q^ten.ra) = / ]T O^i'^e-^ 3 ™-^, (16) 

and take into account the relations between and given by Eq. 1|12|) . The general ansatz for the Fourier 
transformed Q-matrices, 

Q#(k,wi,«a) = /3(27r) d 5 d (k)^[/3g ab ^ 1 , ^ 2 ,o + <5 afc <Wc 2 ,o#K)] + Q^(k, Wl ,a; 2 ) (17) 
consists of the spatially uniform mean-field part and the fluctuating spatial component, Q ab . In Eq. (|17|l . 

D(w) = -M/k, (18) 
while the off-diagonal elements of q oh constitute the ultrametric Parisi matrix 13 



{s/si)q EA < S < Si, 
<?£M Si < s < 1, 



in which si = 2y\q EA T '/ 'k, and is the Edw ards- Anderson order parameter (q aa = q EA ). 
We substitute then this ansatz into Eq. I|14|) and obtain that a(iui n ) consists of three parts, 

a(iu n ) =a^(iu} n ) + a^(iu n )+a^(iu n ). (19) 

a^'(uj n ) is given by 

(1) 16e 2 4q EA A q 
a^'(iuj n ) = - — -II(«w n ) (20) 



where 



n(io;„) - / dre^ T [(35(t) - 1] , (21) 



o 



In the derivation above, we used the result, (1/n) ^ a b q ab q ab = (4/3)q EA A q , where A q = q EA — q(s)ds = q EA S\/2 
is the broken ergodicity parameter, that vanishes linearly with temperature. Note, h ad w e assumed that the vector 
potential entered in a time-translationally invariant manner, the factor of —1 in Eq. (|21|l would not be present. As 
a result, the conductivity would diverge at w„ = as in a superconductor. In <j( 2 )(iw„) we collect the terms that 
contain D(cj n ): 

J (3) W = ItY,D 2 M - TJ2 D(cj m )D(uj m +cj n ) - 2q EA D(cj n ) j . (22) 

The remaining term, a^(iuj n ) arises from the spatially-dependent part Q ^ (k, w 2 ) of the Q-matrices. Writing 
the expression for the current, Eq. I|15f) . m two parts, 



Ji(x,r) = -^1^ / T^i E k pM + jE fl W) [^(k^r.^e^—^-^ + cc], (23) 

ab J ^ ' ui,u 2 \ ui 2 ' 



J ^)— xE/3E/4^4^i E (k 1 +k 2 )Q^(k 1)Wl) c,3)(Q» b (k 2) . 2)W3 ))% i ^- ki )- x e ^--)^) 

ab a — ± wi,W2,W3 

we observe that the contribution from Ji(x, r) vanishes as a result of integration over d d (x — x'). The remaining part 
leads to the result that 

» <3, <-> - 

-4fc 2 r^(k, Wl ,^ 2 ; W „)^ b (k, Wl , W2 )^ b (k, Wl , W2 +c i ;„)] . (25) 
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In Eq. 



gf{Kuj 1 ,uj 2 )^{Qf^,uj 1 ,uj 2 )Qf{-^-uj ll -uj 2 )) = - J2 ^(k,wi,wa), (26) 



4 

/i,!/=l,2 



is the exact propagator for the fluctuations of the Q-fields. The first term is the diamagnetic contribution, while the 
second is paramagnetic and can be formally represented by the standard bubble diagrams 14 and r^ fc (k, u>i, uj 2 ; uj n ) is 
the corresponding vertex function. 

We discuss first the contribution a^(uj n ). The explicit frequency dependence of this part is given simply by the 
prefactor Tl(iuj n )/uj n . Should a phase stiffness exist, this prefactor would be simply proportional to l/oj n , which when 
analytically continued would yield the standard electroma gnet ic response for the conductivity of a superconductor. 
However, this is not the case here. The integral in Eq. 1)21(1 is simply /3(1 — Su„,o) effectively removing thus the 
divergence at zero frequency, unlike what would be the case had we assumed that the vector potential entered the 
action in a time-translationally invariant manner. Note that such an expression although not analytic at uo n — does 
not violate causality because it is, nonetheless, analytic in either the upper or lower half planes. Hence, the 0(2) 
quantum phase glass has a vanishing stiffness in the limit oj n = 0, which of course is the physically relevant regime 
for the dc conductivity. It is in this limit that explorations of all available minima are possible. 

To sec this result more systematically, we analytically continue Ti(iuj n ) using a Hilbert transformation. The denom- 
inator of Eq. H2(J[) can be analytically continued trivially, iu> n — > uj + if], where rj is a positive infinitesimal. We write 
the numerator as 

U(iuj n ) = / dre^" r n(T), n(r) = 05(t) - 1 = n x (r) - n 2 (r) (27) 
Jo 

Although Hi (r) is not an analytic function, we can construct its analytical continuati on u sing the conformal invariance 
condition, 5(r) = 8{t + (3). Performing the integration over the first term in Eq. J57J), we obtain that IIi(a;) = (3. 
Because n 2 (r) = 1 is an analytic function, we adopt the spectral representation 

_ i r j^n^ (28) 

valid for Bose systems, where n 2 (e) = n 2 (e) + ?TI 2 (e). This representation is most convenient for constructing 
the analytical continuation 1 ^. Once we know n 2 (e) , we can obtain both the real and the imaginary parts for real 
frequencies using the Hilbert transformation, 

If deri'(e) , . 

n 2 (w = - / 29 

7r J e — uo — i0 + 



Solving Eq. with n 2 (r) = 1 yields 



n 2 (e) = 2ttS(€) sinh ^ = 7r/3e<5(e). (30) 



The real part is determined by the principal value 



n 2 H = p/ ^^ = /3 /m = /3 / 1 - = ° 

, e-uj 



uj^Q 



Obtained in this fashion, the real and imaginary parts of n 2 (cj) formally satisfy the Kramers-Kronig relations. How- 
ever, both are not regular functions. Hence, it is more convenient to treat the real and imaginary parts of n 2 (u;) as 
limits of two analytic functions. For example, from the regular function, 

g^)=[i(^L^ +l -P^\ (31) 

whose real and imaginary parts satisfy the Kramers-Kronig relations, we obtain the correct limit for n 2 (cj = 0) = 1 
simply from g(uj = 0) = 1, and for uj ^ the limiting procedure, lim^o g(u>) — n 2 (w ^ 0) = 0. As a result, the 
limits, uj = 0, rj — > and rj = 0, lo — * do not commute, a fact which must be considered when we construct the uj = 
conductivity. The correct order of limits is 7/ — > 0,u> — 0. Nonetheless, the advantage of writing n 2 (w) in this fashion 
is that for any non-zero 77, the real and imaginary parts of this g{uo) obey the Kramers-Kronig relations. Combining 
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this representation with 111 (uj) — 1 and iu) n — > to + iS 7 we obtain the analytically continued form for the frequency 
dependence of the conductivity 



*(t\ + 6) . uj 3 — r/Suj 



(if + UJ 2 )(t] 2 + UJ 2 ) (t] 2 + UJ 2 )(S 2 + UJ 2 ) 



w = 0, lim5_ >0 ,7 ) -.o 
cj 7^ 0, lim5_,o,7)^o 



Recall, the correct w = limit is recovered by setting u = and then taking the limit, ?7 — » 0. We find then that the 
contribution of o^ 1 -* (uj) to the conductivity is purely imaginary. The absence of the real part and, as a result, a fo rma l 
violation of the Kramers-Kronig relations here is tied to the presence of the non-analytic function 5(t) in Eq. I|21|l . 
Such non-analyticity at uj = is permissible because the requirement of causality is analyticity in either the upper or 
lower half planes. 

To evaluate the uj — > limit of (uj) we must analytically continue the difference of the first two terms in Eq. 
(EH- Using Eq. JTHJl we obtainiSJi that 

and is some regular function of the infrared cutoff A w and temperature. We see that the contribution (uj = 0) is 
non-critical and metallic. 

Proceeding to the third term, (uj), we first notice that the exact calculation of the propagator G?£,(k, u)\, U2), 
based on the action Eq. I|1U|) is not possible. However, at the quantum critical point in the Gaussian approximation, 

£vO,^i,cj 2 ) = 12 \ 1 1 | 1 1 — 1 =0o(k,u;i,W2), (33) 
k z + \uji\ + \uj 2 \ 

and hence is indepedent of replica and spatial indices. Substitution of this simple replica-symmetric propagator into 
Eq. Ij25(l leads to the zeroth-order result for a^ 3 \uj) as a result of the replica summation. Because the renormalization 
group equations for the coefficients in the action, Eq. (|10(l . lead to runaway to strong coupling for d < d c = 8, it 
is no t p ossible to analyze the behavior of (uj = 0) for the relevant dimensionalities. However, the structure of 
Eq. lf^5|l allows us to make the conclusion that the superconducting contribution of the type p s S(uj) is not expected. 
This can be proven formally by integrating by parts the diamagnetic term and employing the Ward identity. After 
the analytical continuation uj n — > —iuj, we expand the ensuing expression over uj. We obtain that the zero-frequency 
conductivity obeys the scaling form 

albeit the precise form of the function F(x) and, hence, the corresponding temperature dependence can not be 
determined. 

We have obtained an important result that there is no real contribution to the conductivity proportional to p s 5(uj). 
The vanishing of the stiffness is tied to the nature of the vector potential coupling to the glassy order parameter. The 
vector potential couples in a non-time translationally invariant manner to the spin glass order parameter. If, however, 
the system explores only one of the myriad of configurations in the glassy landscapes stiffness appears in agreement 
with the work of Kotliar et. alA However, certainly within a single configuration, the origin of time is irrelevant. 
But this is not the most general case. Quantum mechanically tunneling to all minima is permitted. In this case, the 
stiffness vanishes in agreement with the result^S on the Heisenberg spin glass that the spin stiffness is a transient and 
hence should vanish once tunneling between all minima is present. This result is robust and expected to hold beyond 
the mean-field theory. 



III. BOSONIC CONDUCTIVITY: REPLICA SYMMETRY BREAKING 



Now we generalize our earlier result for the bosonic conductivity. Such a contribution arises only in the case of 
non-zero mean, Jq =/= 0. In this case an ordered phase exists which in the 0(2) case is a superconductor. Hence, in 
the presence of non-zero mean, a new order parameter 

*2(k,T) = (S?(k,T)) (35) 



which is determined by the expectation value of the rotor spin. On the spin glass side of the phase diagram, the 
bosonic excitations of the superconductor develop a mass, m which is equivalent to the inverse correlation length for 
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phase coherence. In the presence of bosonic excitations, the free energy contains the additional terms, 

AF[*,Q] = ]T (fc 2 + ^ + TO 2)|^(k, w „)| 2 



^-Jd d xJ dndr 2 ^(x,Ti)[* b v (x,T 2 )]*Qf v (x,n,T 2 ) 



dr^2[^(x,r)(^(x,r)r] 4 (36) 



At the Gaussian level, with the mean-field spin glass ansatz (Eq. p7jl). the effective Gaussian propagator for the 
bosonic degrees of freedom has the form: 

•fgauss = (k 2 +^l+V\^n\+m 2 )\^ a {k,LU n )\ 2 

-f3q £ <L n ,o^(k,cj n )[^(k,av)]*. (37) 

a,b,k,cj„. 

As we have pointed out previously, the term proportional to 

g ab m 

F gauss cannot be rewritten as an effective mass 
term because this term explicitly couples tp fields with different replica indices. In the case of replica symmetry, that 
is, q ab = qo for all a and b, we have shown that the resultant conductivity is non-zero and given by, 

a hos{LJ = 0,T^0) = -— (38) 

which smoothly crosses over to a = oo in the superconducting state (m = 0). That the bosonic contribution to the 
conductivity should be non-zero is immediately obvious from the |w| term in the action. This term arises entirely due 
to the glass degrees of freedom that naturally provide for dissipation to generate a metallic state. 

We now generalize this result to include replica symmetry breaking. Application of the Kubo formula in this case 
results in a conductivity 



2(e*f ^ f d 2 k r (o) 



a,b,U m 



that depends entirely on the Gaussian propagator for the ip fields. To evaluate this quantity, we need to invert Ea. (|37|l . 
This calculation is difficult to perform for the general type of RSB. However, it can be readily done using the rules 
developed by Mezard and Parish for inverting an ultrametric matrix having a I-step RSB: 



900 



q s < s c 
qi s c < s < 1 



To apply the inversion formula detailed in the Appendix II of Ref. [12] , it is expedient to make the following definitions: 

g= k 2 + v \L\ + m 2 ' ~ 9= (fc2 + ?? K| + i] m ) 1 ( P + ,K| + TO2 )' s ™ = - 2 + ^ * = -.(*- ») (40) 

Application of the inversion formula 18 results in the diagonal 

G = g + pS^Exg 1 —^ + pqoS^.g 2 (41) 



and the off-diagonal elements 



G(s) = /3g ^„,o + ^.oSig — Sc) (42) 



of the propagator. In this representation of the Parisi matrices on the interval [0, 1], the replica indices are absent. 
Nonetheless, a well-defined formula— 



1 ~ f 1 

-TrAB =ab- dsa(s)b(s) (43) 
n Jo 
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exists for taking the trace of a product of two ultrametric matrices A and B, where a and b are the diagonal elements 
of A and B respectively and a and b are the corresponding off-diagonal elements in the continuous representat ion . 
Let's consider here only a simple case of the weak RSB, <§; to 2 , assuming that s c ~ T. Substitution of Eqs. i|41|) 
and l|42|) into Eq. Q39[l . and expanding over /3£i/to 2 results in the following correction to the static conductivity due 
to the replica symmetry breaking 



RSB _2(e*) 2 ^ l-s c [°° xdx _ 4 e 2 



fc = - V-^-ST /„ = 3 ^ (<?1 ~ 9o)(1 ~ Sc) (44) 



Combining this with Ea. l|38|) . we obtain 



Tot _ 4 e 2 ?7 



"bosons - 3^4 fel + (90 - qi)Sc] (45) 

as our total contribution for the bosonic conductivity. If s c = 1, we recover our previous replica symmetric result. For 
the quantum 0(2) spin glass, however, s c oc T, and hence, the correction with s c vanishes at T — 0. Setting s c = 
requires that q(s) — q±. Hence, replica symmetry breaking adds a simple benign constant to the conductivity which 
smoothly crosses over to the replica symmetric result. 



IV. SUMMARY 



We have considered here two separate questions: 1) does the 0(2) vector spin glass have a non-vanishing phase 
stiffness and 2) what is the role of replica symmetry breaking in the bosonic contribution to the conductivity. If the 
answer to the first question were yes, then the answer to the second would be irrelevant as the overall conductivity 
would be infinite. As we have demonstrated clearly, the spin glass order parameter does not provide a superconducting 
contribution to the conductivity at mean field and beyond. Our calculation of the phase stiffness seems to be the 
first based on a direct coupling of the vector potential to the spin glass order parameter which does not assume 
time translational invariance at the beginning. The physical mechanism underlying the vanishing of the spin stiffness 
appears to be the exploration of all configuration minima as a result of quantum tunneling. In addition, we have 
found that replica symmetry breaking provides a small correction to bosonic conductivity. Hence, the bosonic metallic 
state we have found here is robust and represents a clear example of a metallic state in the presence of disorder in 
two dimensions. 
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